Spin polaron theory for the photoemission spectra of layered cobaltates 
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Recently, strong reduction of the quasiparticle peaks and pronounced incoherent structures have 
been observed in the photoemission spectra of layered cobaltates. Surprisingly, these many-body 
effects are found to increase near the band insulator regime. We explain these unexpected obser- 
vations in terms of a novel spin-polaron model for C0O2 planes which is based on a fact of the 
spin-state quasidegeneracy of Co 3+ ions in oxides. Scattering of the photoholes on spin-state fluc- 
tuations suppresses their coherent motion. The observed "peak-dip-hump" type lineshapes are well 
reproduced by the theory. 
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Strongly correlated behavior of electrons is a common 
property of transition metal oxides. This is because the 
bandwidth in these compounds is relatively small com- 
pared to the intraionic Coulomb repulsion between the 
3d electrons. As a result, the celebrated Mott physics [lj 
forms a basis for understanding of a unique properties 
of oxides, such as the high-T c superconductivity and a 
colossal magnetoresistivity. 

Recently, attention has focused on the layered cobalt 
oxides because they exhibit high thermoelectric power 
Q, i.e., the capability to transform heat energy into 
electricity. These compounds consist of a triangular lat- 
tice Co02_planes, separated either by Na layers like in 
NaCo02 [3(] or by thick layers of rock-salt structure in 
so-called "misfit" cobaltates (see 0, HI and references 
therein). Besides controlling the c-axis transport, the 
Na and rock-salt layers introduce also the charge car- 
riers into the C0O2 planes, such that the valence state 
of Co ions is varied in a wide range from nonmagnetic 
Co 3+ t| g 5 = state (as in NaCo02) towards the mag- 
netic Co 4+ ^2 g 5 = 1/2 configuration (in Na^CoO^ at 
small x). 

As the t% g shell of Co 3+ is full, this limit is naturally 
referred to as a band insulator 0, 0], while Co 4+ 5=1/2 
rich compounds fall into the category of Mott systems be- 
cause of unquenched spins. It has therefore been thought 
that layered cobaltates may provide an interesting oppor- 
tunity to monitor the evolution of electronic states from a 
weakly correlated band-insulator regime to the strongly 
correlated Mott-limit by hole doping of NaCo02. Sur- 
prisingly, a completely opposite trend is found experi- 
mentally. The hallmarks of strong correlations such as 
magnetic order 0, strong magnetic field effects 0, 
etc, are most pronounced closer to the Co 3+ composi- 
tions, while Co 4+ 5=1/2 rich compounds behave as 
a moderately correlated metals 0, [§] . The best thermo- 
electric performance of cobaltates is also realized near the 
doped band insulator regime [l(|, thus unusual correla- 
tions in this state and enhanced thermopower are clearly 



interrelated. As a direct evidence of a complex structure 
of doped holes, the recent angular resolved photoemis- 
sion (ARPES) experiments 



ll( observed lineshapes 
that are typical for strongly correlated systems. Para- 
doxically again, the many-body effects in ARPES are 
enhanced approaching the band insulator limit [B[ . 

This Letter presents a theory resolving this puzzling 
situation in layered cobaltates. We show that holes doped 
into nonmagnetic band insulators NaCo02 and misfits 
are indeed a composite objects with a broad energy- 
momentum distribution of their spectral functions. Be- 
sides a reduction of the quasiparticle peaks, they dis play 
also a dispersive incoherent structure as observed [Ellllj. 
The underlying physics behind this unexpected complex- 
ity is in fact quite simple and based on a unique aspect 



of Co + ions, i.e., their spin-state quasidegeneracy, and 
on special lattice geometry of C0O2 layers. 

In oxides, Co 2+ is always in a high-spin 3/2 state while 
Co 4+ ions usually adopt a low-spin 1/2. Roughly, this 
selection is decided by the Hund coupling favoring high- 
spin values, or by the t2 g -e g crystal field splitting IQDq 
supporting low-spin states. An intermediate case is re- 
alized for the Co 3+ ions where the 5 = 0, 1, 2 states 
strongly compete. This leads to a distinct property Co 3+ 
rich oxides called "spin-state-transition" , which is re- 
sponsible for many anomalies such as thermally or dop- 
ing induced spin-state change in LaCoOs El, 

and a 

spin- blockade effect in HoBaCo205.5 14 1, to mention a 
few manifestations of the Janus- like behavior of Co 3+ . 

New situation encountered in layered cobaltates - and 
this is the key point - is that C0O6 octahedra are edge- 
shared. In this geometry, the charge transfer between Co 
ions occurs along the 90° d-p-d bonds, where the largest 
matrix element for the hole motion is that between the 
neighboring orbitals of t2 g and e g symmetry (see t pro- 
cess in Fig. [1]) . This implies that a doped hole dynam- 
ically creates 5=1 i| fl e^ states of Co 3+ fl5| . It is this 
point where the holes - doped into an initially nonmag- 
netic background - become many-body objects dressed 
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FIG. 1: (a) A conventional f-hopping within the low-spin t2 g 
states, (b) f-hopping process creating an excited Co 3+ 5 = 1 
state. Because of the Hund coupling Jh, this triplet state is 
quasidegenerate to the tf s S — state, i.e. the excitation 
energy Et ~ lQDq — 2Jh is low. (c) Bond directions a, b, 
c in the triangular lattice of Co ions and f-coupled orbitals 
on these bonds, (d) Selfenergy of the T spin excitation, (e) 
Selfenergy of the holes in a selfconsistent Born approximation. 



by virtual spin excitations. Formation of spin polarons 
suppresses the plane- wave like motion of holes. However, 
they gain in the kinetic energy by exploiting an addi- 
tional, t hopping channel. 

The layered structure of cobaltates results in quasi 
two-dimensional electronic states like in cuprates. This 
should also lower the symmetry of the t 2g bands 
(threefold-degenerate in the case of a cubic lattice) that 
accommodate doped holes. Indeed, ARPES data shows 
a very simple Fermi surface derived from a single band 
of aig symmetry [H, [n|. This band can be described by 
the nearest-neighbor Hamiltonian H t — — t Ylija- fjafia> 
where the fermionic operators / represent et lff holes. 
They are subject to a conventional Gutzwiller constraint 
(at most one hole at given site) but this is not much 
relevant at small density of holes, i.e., near the band- 
insulator limit. Hence, holes move freely as in a semicon- 
ductor. 

In contrast, the ti g -e g hopping in Fig. [IJb) represents 
the many-body process as it produces the Co 3+ 5 = 1 
excitation. The model describing this process has been 
derived in Ref. Ilfj and reads as follows: 



(described as hole motion). The process, of course, con- 
serves the total spin. In addition, the direction 7 of the 
bond (ij) selects the e g orbital involved in t process [see 
FigJTJc)], according to E g symmetry relations between T 7 

1 \/3 

operators: % = T 3z 2_ r 2, T a / b — — ^T 3z 2_ r 2 ± — T x 2_ y 2. 

Based on this model, we develop a theory for the pho- 
toemission experiments in cobaltates. It is evident from 
§\§ . that by creating and destroying T excitations as they 
propagate, the holes are strongly renormalized and we 
deal with a spin-polaron problem. This resembles the 
problem of doped Mott insulators like cuprates, however, 
the nature of spin excitations is different here because 
of the nonmagnetic ground state. Instead of magnon- 
like propagating modes as in cuprates, fluctuations of 
the very spin value of Co 3+ ions are the cause of the 
spin-polaron physics in cobaltates 17 1. 

For the calculation of the fermionic self-energies, 
we employ the selfconsistent Born approximation [see 
Fig. [T]Je)] , which has extensively been used in the con- 
text of spin-polarons in cuprates 21|. First, we focus 



on spin excitation spectrum. Since a direct e g -e g hop- 
ping in case of 90°-bonds is not allowed by symmetry, 
the bare T spin excitation is a purely local mode, at the 
energy Et- The coupling to the holes in (Q]) shifts and 
broadens this level. Accounting for this effect perturba- 
tively [see Fig. [TJd)] , we obtain the T Green's function 
X> (iu>) = icj — Et — with 



of 2 
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Here, Q$ is the bare electron propagator ^0(^1 le ) = ( le — 
Cfc) 1 with the a ig dispersion on a triangular lattice = 
— 2t(c a + q, + c c ) + fx, where c 7 = cos/c 7 and fc 7 are the 
projections of k on a,6,c directions. The underlying E g 
symmetry of T operators involved in t hopping results 
in the factor Tk — c? a + c| + c 2 c - c a Cb - CbC c - c c c a . We 
neglected a weak momentum dependence of for the 
sake of simplicity. This is justified as long as Et is small 
compared to the spin gap Et- 

Further, we approximate Tk by its Brillioun-zone aver- 
age 3/2, obtaining the simple expressions for Et in terms 
of bare fermionic density of states Nq(x) — ^2 k 6(x — 

,0 

lmE T (E) = -7rf 2 / dxN (x)N (x + E) , (3) 

J-E 



(1) 



Here, T is the spin-triplet excitation generated by hop- 
ping of an electron from Co| + to the e g level of COj + 



ReE T (£) = -V 



Ax 



dy 



2 N (x)N (x + y) 
y 2 -E 2 



(4) 



These equations determine the renormalized spin- 
excitation spectrum pt(E) = — 7r _1 ImI?(ia; — > E + iS) 
used below for calculation of the fermionic selfenergy. 
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The selfenergy diagram in Fig. QJe) reads as: 
2t 2 

Sfc(ie) = --=- P*'X>(-iw) + T k V(iiu)} Q(k', ie+iu 

k ,iu 



where Q ie) = ie— — Sfe(ie). We can write Sfe(w) 
2t 2 [$(w) + r fe S(w)], where 



(5) 



= / dEp T (E) / da; 
Jo JO 



N(x) 



lu — E — x + i<5 



, (6) 



/ dEp T {E) / dx 

JO J-oo 



AT(a;) 



cj + — a; + i(5 



(7) 



The full local density of states N(E) = £ fc A(k, E) and 
its Eg symmetry part N(E) — J^k ^kA(k, E) are func- 
tions of the selfenergy itself, via the spectral functions 
A(k, E) = — 7r _1 Im(?(ie — > E + id). The above equations 
are thus to be solved self-consistently. 

Next, we would like to test the reliability of approxi- 
mations made. To this end, we performed an exact di- 
agonalization of our model H — H t + on a 
hexagonal cluster [Fig.^d)]. We inject one hole on such 
a cluster and calculate its spectral function 

A(k,E) = -hm(GS\fkAz-H)- 1 fijGS) , (8) 

where z = E + Eqs + iS. Since the groundstate \GS) 
contains no holes and the single hole on the cluster rep- 
resents 1/7 doping, the effective doping is about 10%. 
Periodic boundary conditions allow us to access A(k, E) 
at two non-equivalent k points [see Fig. EJd)] . To obtain 
continuous profiles of A(k,E), we have broadened the 
excited states. 

Our model has two parameters: t/t and Ex/t. In fact, 
the ratio t/t ~ 3 follows from the relations t = 2to/3 and 
t/to = tajt-K ~ 2 (where to = t^t^/Apd) 16}- Thus, we 
set below t/t = 3, leaving Et as a free parameter. 

We find that the above equations give results consis- 
tent with the exact diagonalization data, even at rather 
small spin gap values Et ~ t, as shown in Fig. [2ja,c). 
Both approaches lead to spectral functions with a renor- 
malized quasiparticle (qp-) peak whose spectral weight is 
transferred to a pronounced hump structure. [A peculiar 
momentum dependence of the matrix elements r& (note 
that rfc=o = 0) reduces the effect at k — T point]. Sev- 
eral maxima on the hump reflect the presence of multiple 
triplet excitations created by the hole propagation. All 
these are the typical signatures of polaron physics. Mul- 
tiplet structure of the hump will in reality be smeared 
by phonons which are naturally coupled to the t tran- 
sition involving also the orbital sector. Although ex- 
periments 22( indicate that electron-phonon coupling is 
moderate in cobaltates, it may substantially enhance the 
spin-polaron effects, as in case of cuprates 23 1. 
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FIG. 2: Comparison of results from the diagrammatic cal- 
culation at 10 % doping (left) and from exact diagonalization 
(right), (a) Spectral functions and a quasiparticle weight Z 
obtained from Eqs. ©-(J7J) with bare Et ~ 1.2i (renormal- 
ized to Et /t ~ 1 by interactions) and (b) the corresponding 
density of states N(E). The incoherent structure dominates 
N(E). A peak around 2t corresponds to the van Hove singu- 
larity smeared by interaction effects, (c) Spectral functions 
from exact diagonalization at Er/t — 1. (d) The 7-site clus- 
ter in direct space (the dashed line defines the supercell) and 
in reciprocal space (full line is the Brillouin zone boundary). 
The allowed k — T,X points are indicated by (•) and (o). 
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FIG. 3: (Color online) Intensity map of the spectral density 
A(k, E) of the Co 4+ holes along M-T-K path in the Brillouin 
zone calculated at 30 % doping and Et = It. As the hole en- 
ergy reaches the renormalized spin excitation energy Et, the 
qp-peak broadens and its weight is transferred to a broad, in- 
coherent structure. This results in a "peak-dip-hump" profile 
of A(k, E) seen also in Fig. [2] The top of the smoothed hump 
structure is indicated by triangles. The dashed line shows the 
bare dispersion. The Fermi surface is shown in the inset. 
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FIG. 4: Top panel: Full energy width at half maximum of 
the quasiparticle peak along the M~r dispersion curve (see 
bottom-right panel) plotted as a function of the binding en- 
ergy. The part from /cf to M is accessible by ARPES exper- 
iments. Strong qp-damping below ~ — 1.4£ is due to a scat- 
tering on S = 1 excitations. The spectral weight of the qp- 
peak obtained by a direct integration is indicated by squares. 
When damping is small, it coincides with a conventional qp- 
residue (1 — dT,' / duj)" 1 . Lower panel: The T-exciton spectral 
function (left), and renormalized hole dispersion (right). 



I At for the qp-damping (Fig. [4| is then w 140 meV, in 
nice agreement with experiment (see Fig. 2c of Ref. . 
Physically, dilute spin-polarons are expected to be pinned 
by disorder, thus qp-peaks should be suppressed at low 
hole doping. 

To summarize, we have presented a theory for the pho- 
toemission experiments in layered cobaltates. A strong 
damping of quasiparticles, their reduced spectral weights, 
broad and dispersive incoherent structures, and a "peak- 
dip-hump" type lineshapes all find a coherent explana- 
tion within our model. We thus conclude that unusual 
correlations observed near the band insulator regime of 
cobaltates are the direct manifestation of the spin-state 
quasidegeneracy of Co ions, and this intrinsic feature 
of cobaltates should be the key for understanding of 
their unique physical properties, such as the high ther- 
mopower. In particular, the spin-polaronic nature of 
charge carriers should be essential for explanation of its 
remarkable magnetic field sensitivity 0]. 

We would like to thank B. Kcimer for stimulating dis- 
cussions. This work was partially supported by the Min- 
istry of Education of CR (MSM0021622410). 



To illustrate the gross features of the hole renormal- 
ization, in Fig. [3] we show a complete map of the spectral 
function along M-F-K path in the Brillouin zone. We 
have used a representative value Et = 2< which is renor- 
malized by holes to Et ~ lAt. Compared to the bare 
dispersion, the bandwidth of the renormalized holes is 
reduced by a factor ~ 2. The main observation here is 
that as the hole energy reaches Et, the dynamical gen- 
eration of S = 1 excitations becomes very intense and a 
broad incoherent response develops, leading to the pro- 
nounced "peak-dip- hump" structure of A(k, E). Follow- 
ing the maximum of the smoothed hump structure, we 
observe its strong dispersion (stemming also from inco- 
herent t hopping). 

Fig. [3] suggests a possible determination of Et from 
the quasiparticle damping. To address this problem, in 
Fig. Q] we show the energy width of the qp-peak following 
its dispersion curve. The sharp onset of the damping 
at the binding energy ss —lAt is clearly related to the 
maximum of the spin-excitation spectral function pt{E). 
In addition, Fig.[4]shows the weight of the qp-peak which 
is fc-dependent (mainly due to the matrix element F&). 

Comparison of Figs. [3] and [4] with the data of Refs. 5, 
111 reveals a remarkable correspondence between theory 
and experiment. In particular, both the qp-peak and the 
hump dispersions (see Figs. 2 and 3 of Ref. lfj|) are well 
reproduced by theory, considering t sa lOOmeV suggested 
by the band structure fit [24J. The onset energy Et ~ 
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